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The theory of epitaxial strain energy is extended beyond the harmonic approximation to account 
for large film/substrate lattice mismatch. We find that for fee noble metals (i) directions (001) 
and (111) soften under tensile biaxial strain (unlike zincblende semiconductors) while (ii) (110) and 
(201) soften under compressive biaxial strain. Consequently, (iii) upon sufficient compression (201) 
becomes the softest direction (lowest elastic energy), but (iv) (110) is the hardest direction for large 
tensile strain, (v) The dramatic softening of (001) in fee noble metals upon biaxial tensile strain is 
caused by small fcc/bcc energy differences for these materials. These results can be used in selecting 
the substrate orientation for effective epitaxial growth of pure elements and A p B q superlattices, as 
well as to explain the shapes of coherent precipitates in phase separating alloys. 

PACS numbers: 62.20.Dc, 68.60.-p, 81.10.Aj 



When a material is compressed hydro statically, its en- 
ergy rises steeply because all three crystal axes are de- 
formed (dashed line in Fig. |l|). When the same mate- 
rial is confined coherently onto a substrate ("coherent 
epitaxy 1 ' 1 ) with lattice constant a s , the energy rises less 
steeply (solid line in Fig. |l|) since it is deformed only 
along the crystal axes in the substrate plane and allowed 
to relax (and thus, lower its energy) in the third direction 
G. This "epitaxial softening" can be quantified by the 
dimensionless parameter 
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giving the ratio between the epitaxial increase in energy 
due to biaxial deformation to a s , and the hydrostatic in- 
crease in energy due to triaxial deformation to the same 
a s . Because the biaxial strain energy AE epl (a s ,G) de- 
pends on the strain direction G, so does q(a s ,G). In 
growing coherent epitaxial filmsjj it is desirable to min- 
imize AE cpl (a s ,G) [or, equivalently, for a fixed a s mini- 
mize q(a s ,G)], so as to avoid or reduce dislocations and 
other strain-induced film/substrate defects. It is hence 
important to select substrates a s and growth directions G 
that entail miniu^Lslxainjjpncrgy. Harmonic continuum 
elasticity theoryaBuoQuffl makes definitive predictions 
for the a s - and G-dependence of epitaxial strain energy: 
(i) q(a s ,G) does not depend on a s , and (ii) the "softest 
direction" is (001) if A = [C 44 - |(Cn - C 12 )] > 0, while 
if A < then the softest direction is (111). For most 
fee metals and semiconductors A > 0, whereas A < 
for ionic salts (PbS, AgBr, NaCl, KC1) and several bee 
metals (Nb, V, Mo, Cr). The selection of substrate ori- 
entation for many years has been guided by these predic- 
tions of harmonic elasticity, summarized compactly by 
the expressionQiJ 
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where B = |(Cn + 2C12) is the bulk modulus and 

/harm 

(G) is a purely geometric function of the spherical 
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FIG. 1. Epitaxial (AE cpl (a s , G), solid line) and hydrostatic 
(AE hulk (a s ), dashed line) deformation energies for fee Cu. 
The dotted line shows a fit to AE cpl (a 3 ,G) using 3-rd order 
polynomial in a s . 

angles formed by G: 
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Ki(4>,0) are Kubic harmonicstLHl (linear combinations of 
spherical harmonics forming irreducible representations 
of the cubic group), and the spherical angles </> and 9 are 
measured with respect to Cartesian axes oriented along 
the edges of the conventional fee cubic cell. Figures ||(b) 
and @(c) show (7harm(G) of Au and Cu, demonstrating 
that (001) and (111) are, respectively, the softest and 
hardest directions. 

One expects to find corrections to the harmonic be- 
havior predicted by Eq. (||) due to the fact that strong 
deformations may alter the electronic charge density (and 
thus, the elastic response) beyond the validity of conven- 
tional harmonic elasticity theory. However, the theory of 
elasticity itself does not indicate when and how expres- 
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(a) Au, a s =3.56A 
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FIG. 2. Parametric plot of the epitaxial softening func- 
tion q(a a ,G) vs. substrate orientation G for (a)-(b) Au and 
(c)-(d) Cu, at a few values of the substrate lattice constant 
a s . The equilibrium (eq) lattice constants of fee Au and fee 
Cu are shown in (b) and (c), respectively. 

sion (||) will fail for a given system. Figure [3] shows that 
at least a fourth-order polynomial in a s is needed to re- 
produce the qualitative structure in AE cpl (a s , [001]) for 
Cu, as obtained directly from accurate electronic struc- 
ture calculations (see below). The three extremal points 
in AE epl (a s , [001]) suggest that a simple extension of 
the harmonic elasticity to include third-order elastic con- 
stants is not sufficient (the dotted line in Fig. [I] shows the 
behavior of the strain energy predicted by a 3-rd order 
polynomial fit to a few data points around a cq = 3.56 A). 
The goal of the present work is to propose a simple an- 
alytic formula for q(a s ,G) [and therefore AE epl (a s , G)] 
which works for arbitrary strain magnitude and direction 
G. 

We have calculated AE epl (a s ,G) and q(a s ,G) for Ag, 
Au, Cu and Ni along six principal directions (001), (111), 
(110), (113), (201) and (221) using the local density ap- 
proximation (LDA) JLH as implemented rW the linear aug- 
mented plane wave (LAPW) methodJla This approach 
gives the total energy as a function of biaxial distortion, 
including all anharmonic effects. We find that our results 
can be fitted accurately to a generalization of the Kubic 
harmonic expansion of Eq. (|^) : 
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where bi(a s ) is now a function of a s , and the angular 
momenta extended beyond the harmonic limit (/ = 0, 4) 
to include I = 6, 8, 10, . . . The analytic representation of 
Eq. (||) allows us to explore the elastically soft and elas- 
tically hard directions for any strain, even outside the 
regime of harmonic elasticity. We find that for fee noble 
metals q(a s , [001]) and q(a s , [111]) soften as a s expands 
(tensile strain) while q(a B , [110]) and q(a s , [201]) soften 
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FIG. 3. The calculated epitaxial softening functions 
q(a 3 ,G) vs. substrate lattice parameter a s for Ag, Au, Cu 
and Ni. Points represent the directly calculated LDA values 
and lines show the fit using Eq. (m. 



as a s contracts (compressive strain). Consequently, upon 
sufficient compression (001) is no longer the softest direc- 
tion, but rather (201) is (e.g., Ag, Au). Also, upon suffi- 
cient expansion, (111) is no longer the hardest direction, 
which is now (110) (e.g., Cu, Ni). Finally, the dramatic 
softening of q(a s , [001]) upon biaxial expansion reflects 
the existence of low-energy "excited" bec and bet crys- 
tal forms. These results can better guide the selection of 
substrate orientation for effective epitaxial growth. 

We have performed LDA calculations for biaxial com- 
pression (a s < a eq ) of Au and Ag, for biaxial expansion 
(a s > a cq ) of Ni, while for Cu we considered both expan- 
sion and compression. Table [j] gives the calculated LDA 
elastic constants, which are typically 10-20% higher than 
the experimental values for Cu, Ag and Au. Larger dis- 
crepancy in the case of fee Ni is due to the neglect of spin- 
polarization effects in our LAPW calculations. Figure || 
shows a line plot of q(a s , G) vs. a s for a few directions 
G, while Fig. ^| shows a parametric plot of q(a s ,G) vs. 
G for a few substrate lattice constants a s . We find sev- 
eral surprising predictions relative to the expectations of 
harmonic elasticity. Although the harmonic expression 
9harm(G) [Eq. ^ is exactly satisfied at the equilibrium 



lattice constant a s = a cq [where lim as _ >aoq bi(a s ) = in 
Eq. (Q)], this ceases to be so as a s deviates from a cq . 
This is apparent from the dependence of q on a s (absent 
in the harmonic theory), from the crossing of q for dif- 
ferent G values, and from the development of new lobes 
and minima in Fig. |^ with the change in a s . Such effects 
occur at Aa/a < 4%, suggesting a rather small range 
of validity of the harmonic approximation. Furthermore, 
in the harmonic elasticity theory, if (001) is the softest 
direction then (111) must be the hardest direction, and 
vice versa. Figure |3| shows that this does not hold for suf- 
ficiently deformed films: the hardest direction in Ni and 
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TABLE I. The calculated and experimental (Ref. |l7[ ) elas- 
tic constants of fee Ag, fee Au, fee Cu and fee Ni (in Mbar). 
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Calc. Expt 




1.52 


1.24 


2.14 1.86 


2.30 1.68 


3.30 2.47 


Cl2 


1.09 


0.93 


1.73 1.57 


1.58 1.21 


2.20 1.47 


C44 


0.61 


0.46 


0.37 0.42 


0.99 0.75 


1.36 1.25 



Cu for a s ^> a cq is (201), while the hardest directions for 
Ag and Au at a s <C a oq are (111) and (001). 

The new results for fee noble metals, apparent from 
our self-consistent LDA calculations, are: 

(i) q(a s , [001]) and, to a lesser extent, q(a s , [111]) soften 
as a s expands (tensile biaxial strain). 

(ii) q{a s , [110]) and q(a s , [201]) soften as a s is com- 
pressed. 

(iii) As a result of (i) and (ii) above, we find that 
upon sufficient compression, (001) is no longer the elasti- 
cally softest direction, but (201) and, to a lesser extent, 
(110) are. The hardest direction upon compression is still 

HI • 

(iv) Upon sufficient expansion, (111) is no longer the 
hardest direction, but (110) is (Cu, Ni). The softest di- 
rection upon expansion is still (001). 

We find that result (i) is a reflection of the existence of 
a low-energy bec and bet "excited" structures. IndeetL 
(001) strain applied to fee lattice defines a Bain path,E3 
transforming fee into bec via the body-centered tetrago- 
nal structure. This intermediate structure is character- 
ized by tetragonal unit cell dimensions a and c. When 
minimizing the total energy of a (001) biaxially deformed 
solid with respect to the out-of-plane lattice vector c at 
each in-plane lattice parameter a s , one finds an "epitax- 
ial (001) Bain pattt'Eil (solid line in Fig. p]). For noble 
metals having the fee structure (c/a = V2) at equilib- 
rium, this deformation path contains the bec (c/a = 1) 
saddle point, audpthci.bct (c/a = 0.96 in the case of Cu) 
local minimum. oEjll3 The low amplitude of the epitaxial 
(001) Bain path relative to the hydrostatic path (Fig. [l]) 
defines the softness of q(a s , [001]) via Eq. ([jj), and there- 
fore is a direct manifestation of the small fcc/bcc and 
fec/bet energy differences. Indeed, the epitaxial soften- 
ing function at c/a = 1 is given by 
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where A£; bcc / fcc = Et 



= (V c b q cc ) - Et t (V^ c ) 
noble metals V/? q c « Ksq> the bec point is reached at 
a s 2^a oq (see Fig. |). 

A similar argument explains the softening of 
q(a s , [111]), since the epitaxial (111) Bain path also con- 
nects cubijU-.fcc (c/a = -y/6) and bec (c/a = -s/6/4) 
structures.Ej This property of the (001) and (111) Bain 
paths is caused by the cubic symmetry of the crystal 
at the bec point c/a = 1, requiring it to be a local 
extremum of the total energy. It is interesting that 
in zincblende solids (two interpenetrating fee lattices), 
tetragonal (001) expansion does not lead to an extremal 
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point in AE cpl (a s , [001]) since the crystal at the "bec 
point" c/a = 1 does not posses higher symmetry than 
at c/a 5^1. As a consequence, AE cpl (a s , [001]) is a 
monotonically increasing function of a s , and (001) tensile 
strain does not lead to a softening of q(a s , [001]). 

The consequences of our findings are as follows: (a) 
Films of fee noble metals under tensile biaxial strain pos- 
sess the lowest strain energy for strain direction (001), 
while large compressive biaxial strains will have lower 
strain energy for (201). (b) When growing an S p L q su- 
perlattice with components having small (S) and large 
(L) lattice constants (e.g., S = Cu, L = Au), the system 
rich in L (q p) will have a low (001) elastic strain en- 
ergy due to easy expansion of S, while a system rich in S 
(q -C p) will have a low (201) elastic energy due to easy 
compression of L. 
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